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We introduce a Cellular Dynamical Cluster Approximation (CDCA) to study the nature of the 
Mott insulator-metal transition in the extended Hubbard model on a square lattice. At strong 
coupling, a d-wave Mott insulator is obtained. Hole doping drives a first order Mott transition 
to a non-Fermi (nFL) liquid metal. Remarkably, this nFL is caused by an Anderson orthogonality 
catastrophe at low energies due to the non-trivial competition between strong, non-local interactions 
and hopping. This constitutes the first explicit realisation of Anderson's Luttinger liquid idea in 
two dimensions. Many experimental responses in the "strange metal" phase found around optimal 
doping in cuprates are understood naturally within our approach. 

PACS numbers: 71.30.+h,72.10.-d,74.20.Mn 



The problem of high-T c superconductivity (HTSC) in 
quasi-two-dimensional cuprates remains an outstanding 
problem of modern condensed matter physics. HTSC 
arises from a new, non-Fermi liquid (nFL) state, charac- 
terised by highly anomalous normal state response of the 
underdoped cuprates, continuing up to optimal doping, 
x pt (where the SC T c is maximum). Beyond x op t, a low- 
T Fermi liquid state seems to be recovered Q. (i) Super- 
conductivity itself is unconventional, and of d-wave sym- 
metry [2j. Subsequently, it was also found that the nor- 
mal state pseudogap (in the one-particle spectral func- 
tion) possesses d-wave symmetry. Given the problems 
encountered in relatingthe PG phase to precursor pairing 
fluctuations above T c [3|], a new broken symmetry state, 
with d-density wave (d-DW) order, was proposed the 
d-DW state was hypothesised to compete with d-SC, be- 
ing favoured in the underdoped region above T c , (ii) 
On the other hand, careful measurements [l[ reveal that 
the cuprates remain insulating (at very low T if one sup- 
presses d-SC by high magnetic fields) almost up to opti- 
mal doping. Given the amount of disorder involved, this 
insulator would be a Mott-Hubbard- Anderson insulator. 
If the Mott transition occurs around x op t , a natural ques- 
tion would be to inquire about its relation to the maxi- 
mum of T c {x) at the same point. Does the Mott insula- 
tor then have d-wave order beginning from the undoped 
(x = 0) state? (AF/SC long range order is precluded in 
2D by the Mermin- Wagner theorem, but d-wave order is 
not). How does this evolve to a nFL metallic state upon 
hole doping? 

An attempt to reconcile both (i) and (ii) above would 
require that the Mott transition is an insulator-metal 
(IM) transition from a d-wave insulator to a d-wave 
metal. This gives rise to the question: Is the quantum 
critical point (QCP) invoked to understand cuprates [f| 



the one associated with such a. T — Mott transition? 
(Recall that the T — x phase diagram shows that AF-Neel 
order is destroyed well before T c (x) achieves its maximum 
around x op t)- If so > one would expect strongly diver- 
gent, d-wave charge/spin fluctuations around this tran- 
sition (around x op t) to naturally mediate high-T c SC of 
d-wave type. This would constitute a very new picture 
of the physical response of cuprates. 

Additionally, strong dynamical spectral weight trans- 
fer (SWT), a ubiquitous feature of correlated systems, 
is exhibited clearly [j| as the Mott insulator is doped 
(this occurs over a region of 0(4.0 eV), and can only 
be caused by electronic correlations). Theoretically, non- 
local effects giving rise to d-wave correlations can be ac- 
cessed only by cluster-DMFT (see review in [6]). Such 
studies indeed find good qualitative agreement with the 
T — x phase diagram for the cuprates. However, to our 
best knowledge, the compelling link between the unusual 
nFL responses, the T = Mott transition around x op t, 
and the high-Tc SC of d-wave type still remains to be 
described consistently. 

Here, with this motivation, we study the Mott transi- 
tion in a two-dimensional model (without chemical disor- 
der) for cuprates. We construct a new cellular-dynamical 
mean field approximation (C-DMFA) and use it to study 
the doping-driven Mott transition from a d-wave Mott 
insulator to a d-wave non-FL metal. We consider the 
extended Hubbard model (EHM) in the strong coupling 
limit as an effective model for Cu-0 layers in cuprates 0] • 
The Hamiltonian is 

H = -t ( c L c j° + h.c) + U ^n iT n a 
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where n, = and J7 >> V, t, and -Ef is the chem- 

ical potential. 

Next, we formulate a new Dynamical Cluster Approx- 
imation, which wc christen Cellular Dynamical Cluster 
Approximation (CDCA). Given indications from several 
experiments @, our aim is to treat inter-site elec- 
tronic correlations exactly on a single plaquette (2x2) 
and consider the inter-plaquette correlations (resulting 
from carrier hopping and n.n interactions) within the 
CDCA. This is done by mapping the original EHM to 
an electronic model describing coupled plaquettes, and 
is a multi-stage procedure as described below: (1) Con- 
sider the EHM on a single 2x2 square plaquette, C. In 
the limit U — ► oo, the EHM reads 



Hubbard model with U — ► oo, and a fictitious "crystal 
field" term which raises the d-wave state above the dou- 
bly degenerate p XtV and extended-s states. (3) We now 
connect C to other plaquettes via the hopping, i, and n.n 
coupling, V, in the original EHM. To do this, we split the 
2D lattice into plaquettes covering the whole lattice. The 
hopping and n.n coupling between neighboring sites on 
different neighboring plaquettes (in terms of original site 
fermions) are then expressed in terms of the plaquette 
fcrmions (a^u). In the 2D square geometry, V,t acting 
between n.n plaquettes means that sites (2, 3) on plaque- 
tte p, are connected to sites (1,4) on (/i + 1) along x, 
while (1, 2) on fi are connected to (3, 4) on (// + 1) along 
y. Namely, the Hamiltonian coupling the given plaquette 
fi to the rest of the lattice is 



H = —t J2 (Xlx ja + h.c) + v J2 n ^ 

<7,<i,j> = l <i,j> = l 

4 

- E F n io , (2) 

<T 7 i— 1 

where the no-double occupancy constraint resulting from 
the U — ► oo limit is implied, and appears as a local con- 
straint on H. Using the Gutzwiller approximation Q, 
this becomes an unconstrained fermion model, but with 
V = V(x),t — t(x). (2) Construct the four possi- 
ble symmetry-adapted linear combinations of the site 
fermions to "plaquette-centered" fermions: a^ a — Ui^Ciu 
with i — [1,2,3,4] for a 4-site square. Explicitly, a\„ = 
i(ci cr +C2 CT +C3 CT +C4 Cr ), with a 2(T = ^{c la -c 2r j+c 3a -C4 Cr ), 
and a 3o . = \{c la + c 2a - c 3tx - c ia ), and finally, a 4(T = 
\{c\ a — C2a — C3a + C4 a ). Notice that ai transforms 
as an s-wave operator, a 2 as a d-wave operator, while 
0,3,0,4 transform like p-wave operators. The one plaque- 
tte Hamiltonian becomes (fi is centered on C) 



He = -2t^(n a i MCT - n a 2 M cr) + V 53 n a ^ a n afJ _p 

— Ep 2_j n at*a , (3) 



Thus, for half- 



where a, (3 = [1, 2, 3, 4], and 
filling, t < 0, and U — » 00, the exact spectrum of the 
EHM on an isolated square consists of four singly occu- 
pied plaquette states, E\ = —It, E%$ — 0, E4 = 2t, which 
we identify with the "lower Hubbard band" states, and 
the corresponding states raised upwards by V , which we 
identify with the "upper Hubbard band" (charge trans- 
fer) states. In analogy with the usual Hubbard model, the 
UHB states correspond to high-energy states describing 
the dynamics of an electron hopping onto a given local 
"site" (plaquette) in a specified "orbital" when the other 
"orbitals" are already occupied. He in the plaquette ba- 
sis is exactly similar to the local limit of a four-orbital 



H' t = -*^( C W C2 ^+ e «"°- + C t^ C 3,fi+e x ,<T + h.c) 

- ^(iAc+^^ + ^^rt^+M. ( 4 ) 



and 



H' V = U^(n 2 , M «4, M +e* + "3,M n 4, M + e J 

+ V^2(n 4 ^n hfl+ey + n 2 , IM n 3tlM+ey ) . (5) 
In terms of the plaquette fermions, H[ + H' v = H' 



becomes 



H' = - 53 tap{a} taa aii,+e tt , v i3<, + h.c) + — 53 n m n^ + x, 



after 



16 



H,af3 



V 



— l t t 

g 2-^1 ^ a lp.<y a 2fj,a-' a 2,p.+e X} ya-' a l,fJ.+e :Ciy a 

+ 1, 2 — ► 3, 4 + h.c) + other terms 



(6) 



where v) now denote neighboring clusters. Here, 
t a , a = (i/2)(l, -1,1,-1) for a = [1,2,3,4], and i 13 , 2 4 = 
t/2(l, — 1) along x and t/2(— 1,1) along y. Our to- 
tal Hamiltonian in terms of plaquette centered fermions 
is multi-orbital Hubbard model with the intra-orbital 
Hubbard U — > 00, the "inter-orbital", plaquette-local 
Coulomb interaction V, along with an explicit non-local, 
interaction term. This term connects neighboring clus- 
ters, i.e, it is of 0(1/D 2 ). We decouple it in the gen- 
eralised Hartree-Fock scheme (exact to 0(1/ D 2 )). The 
other terms are unimportant near half-filling (they couple 
ex-s and p XtV states), and are neglected henceforth. 

Two d-wave order parameters are obtained: the 
particle-hole average, A dl = (a| Mcr a 2Mcr + a-l^a^a) and 

the d-wave pairing average, A12 = { a \ tia a \iy-a) ■> as re ~ 
quired. A multi-orbital dynamical mean field approxima- 
tion (MO-DMFA) for this Hamiltonian is exactly equiv- 
alent to a cluster-DFMA for the original EHM on the 2D 
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FIG. 1: The "orbital resolved" and the total many-particle 
density of states (DOS) for the 2D EHM at half-filling, 
(n) = 1. Notice the appearance of "orbital" dependent Mott- 
Hubbard gaps, and the d-wave insulator, as explained in the 
text. 



FIG. 2: Same as in Fig.(l), but for hole doping, x — 0.1 
(dashed lines) and x = 0.2 (solid lines). The "orbital selec- 
tive" character of the Mott transition is explicitly seen, as is 
the fact that the ex-s-wave DOS remains Mott localised in 
the metallic phase. 



square lattice. Thus, we solve H = He + H' using the 
MO-DMFA. 

In order to make progress, we use insights gained from 
detailed analyses of the general multi-orbital problem in 
usual DMFA [9. llOfl. Here, multi-orbital models without 
the inter-orbital, intersite pairing terms were solved us- 
ing IPT and QMC, both leading to qualitatively similar 
conclusions. The Mott insulator has ferro-orbital (in our 
case d-wave) order: this is a d-wave Mott insulator. Away 
from n = 1, the highest-lying (d-wave) band becomes 
metallic, while the lower-lying (here, one or more of ex- 
s,Px,y) bands remain Mott insulating. For our model, 
this would describe a first order Mott transition from a 
d-wave insulator to a d-wave metal for any x = (1— n) > 
in absence of chemical (doping-induced) disorder. As we 
show below, this ci-wave metal is not a Fermi liquid. 

We solve our cluster (mapped 4-orbital) problem using 
multi-orbital iterated perturbation theory (MO-IPT). As 
discussed in recent [llj work, the MO-IPT turns out to 
be an accurate solver for multi-orbital problems at ar- 
bitrary fillin gs, n < 1. We refer the reader to relevant 
references (il llll] for methodical details. We focus on the 
Mott-Hubbard I-M transition as a function of hole dop- 
ing. The role of the d-wave p-h (PG) and SC orders are 
described in an accompanying work [181 ] . Here, we focus 
on the region of the T — x phase diagram for cuprates 
without d-wave PG/SC order; i.e, on the strange metal 
phase terminating at T — at x opt - 

We now describe our results. In Fig. [IJ we show the 
"orbital" resolved and total one-particle spectral func- 
tion for our model with V = 4.0 eV, W = 8t = 4.0 eV 
for half-filling. First, it is clear that the occupied states 



closest to Ep(= 0) have d-wave symetry (observe that 
d ® s' — d,p x ® p y = d). The Mott-Hubbard gap, 
Ajuff * 1.0 eV. Very interestingly, different "orbital" 
DOS have different Mott gaps: this is a benchmark fea- 
ture of the (d-wave) ground state breaking the lattice 
point-group symmetry. In terms of original c-fermions, 
these features result directly from an explicit (and ex- 
act to 0(1/1?)) consideration of the competition between 
strong, intersite correlations and one-electron hopping in 
the strong coupling (U >> V,t) limit of the EHM. To- 
gether with the rich Hubbard band structures in Fig.[TJ it 
describes the internal structure of the Mott-Hubbard gap 
(inaccessible in d = oo). We predict that polarised op- 
tics, XAS, Raman measurements would show anisotropic 
spectral features reflecting those found in the DOS: this 
may already have been seen. 

To extend our calculation to n < 1, the MO-DMFT 
equations need to be constrained by the Friedel-Luttinger 
sum rule, An = 4(1 — a;) = J^Z, ptotai(u)du> for our cluster 
problem. Ep is then determined self-consistently within 
the MO-DMFT loop at each stage using ptotai{u) com- 
puted in each iteration, and the process is iterated to 
convergence. The spectral functions for n < 1 are shown 
in Fig. [2 As expected generally, hole doping trans- 
fers dynamical spectral weight over large energy scales 
O(4.0 eV) from the Hubbard bands into the Mott gap. 
The MIT is discontinuous, with the carrier density, (n), 
jumping discontinuously from zero to a finite value off 
half- filling (not shown), as generically expected within 
MO-IPT [9(. What is more interesting, and crucially 
important, is that this SWT is also strongly "orbital" 
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dependent: clearly, our results show that the highest- 
lying d-wave DOS becomes the most metallic, the p XtV 
DOS much less so, while the ex-s-wave DOS remains 
Mott insulating. In "orbital" language, this is an "or- 
bital selective" Mott transition, discussed in detail in 
recent LDA+DMFT work on real oxides @, Q. In our 
CDCA applied to the 2D EHM, this feature directly fol- 
lows from the strong, dynamical scattering between elec- 
tronic channels with different "orbital" symmetries in a 
situation where these are split on the single cluster (rem- 
iniscent of crystal-field splitting in usual multi-orbital 
models [I [HI). 

A very remarkable fact now emerges. Since the ex-s- 
DOS remains Mott insulating in the metal, strong scat- 
tering between the quasi-itinerant <i-DOS and the lo- 
calised ex-s-DOS results in an exact mapping of our 
problem to an effective Falicov-Kimball model [^, [Io| . 
with resulting X-ray edge (XRE) behavior. The d- 
wave DOS then shows a power-law fall off at low en- 
ergy in contrast to the smeared lorentzian of a FL metal: 
p d (u) ~ \tu\-^- a \ with a = tan- 1 {V(x)/W(x)). Hence, 
so does the total DOS. The corresponding self-energy 
is S(cj) ~ \to\^- a \ with a < 1 for V > 12]. This 
implies no FL quasiparticle residue: Zfl = 0, and the 
metallic state is a non-Fermi liquid. Amazingly, this is 
exactly in agreement with Anderson's proposal [131 ] for 
the breakdown of FL theory in the 2D Hubbard model. 
It is indeed remarkable that a close numerical estima- 
tion of the DOS near E F (= 0) gives a ~ 0.2,0.1 for 
x = 0.1,0.2. To our best knowledge, this is the first 
explicit demonstration of Anderson's conjecture within 
C-DMFT schemes. Recall that XRE behavior was sug- 
gested in the past 14| in mult i- channel, single-impurity 
Anderson (and their D = oo lattice counterparts) mod- 
els. 

Armed with this insight, several unusual features of 
optimally doped cuprates are understood naturally. The 
photoemission intensity, as measured in PES, measures 
the one-particle Green's function in the sudden approx- 
imation: hence, IpEg(w) ~ w^ 1- ") as u — > Ep, 
as observed 1131 . The in-plane optical conductivity, 



<7 xx (oj) ~ uj-i 1 - 2 ") from a direct scaling argument, and 
drops like u~°- 6 (x = 0.1), m~°- 8 (x = 0.2) very close to 
that measured experimentally [151 ]. Using an "extended 
Drude" fit [Hj ]. we also obtain that the scattering rate, 
r _1 (oj) ~ ti)' 1_2a ', as indeed observed. Further, the op- 
tical phase angle, </> = tan~ 1 (a2/&i) = (1 — 2a)7r/2, is 
independent of u>, and increases with decreasing x , as 
seen. Finally, the electronic Raman scattering intensity 
in the B\ g channel follows from the local limit of the 
Shastry-Shraiman relation (valid in our MO-DMFT) for 
intraband transitions [3]: Ir(oj) = oj<ji(oj) ~ L0 2a , and 
is very weakly w-dependent, again as seemingly observed. 
Very interestingly, this non-FL behavior extends up to 



rather high energies O(1.0 eV), fully consistent with that 
used phenomenologically in optical analysis [15j . Follow- 
ing Anderson [13], the tunnelling spectrum will also ex- 
hibit power-law anomalies, reflecting that of the DOS. 

What is the dominant instability of this singular 
non-FL phase? Using the XRE mapping, it is clear 
that the "excitonic" correlation function, Xi2"(w) = 
J^dre^ial^ir^al^a^iO)) ~ M"^"" 2 ) is also 
singular. This implies soft d-wave p — h modes, and so 
the dominant instability will be to a d-wave ph pseudo- 
gapped state. Details will be presented elsewhere flij . 

In conclusion, we have formulated a new cluster- 
DMFT to treat the effects of intersite electronic corre- 
lations (exactly) on the Mott transition in the 2D EHM. 
The d-wave Mott-Hubbard insulator found for n = 1 
undergoes a first-order Mott-Hubbard transition to a d- 
wave, non-FL metal upon hole doping. Remarkably, in 
absence of the d-pseudogap, this nFL metal shows low en- 
ergy singular behavior of the one-particle DOS/e?-wave 
"excitonic" correlator, exactly as postulated by Ander- 
son. Our results provide an appealing, microscopic ex- 
planation for a host of normal state anomalies character- 
ising the "strange metal" phase of cuprates near optimal 
doping. 
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